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ABSTRACT 

This report describes the life analysis used to compare copper tubes and milled copper channels 
for rocket engine cooling passages. Copper tubes were chosen as a possible replacement for the 
existing milled copper channel configuration because (1) they offer increased surface area for addi- 
tional enthalpy extraction; (2) they have ideal pressure vessel characteristics; and (3) the shape 
of the tube is believed to allow free expansion, thus accommodating the strain resulting from thermal 
expansion. The analysis was a two-dimensional elastic-plastic comparison, using a finite element 
method, to illustrate that, under the same thermal and pressure loading, the compliant shape of the 
tube increases the life of the chamber. The analysis indicates that for a hot-gas-side-wall tempera- 
ture of 1000 °F the critical strain decreases from 1.25 percent in the channel to 0.94 percent in the 
tube. Since the life of rocket thrust chambers is most often limited by cyclic strain or strain 

range, this decrease corresponds to an expected tube life which is nearly twice the channel life. 


INTRODUCTION 

The future cryogenic Space Transfer Vehicle Engine is required to have a long life and be reli- 
able. This can be achieved with a thrust chamber that is compliant. Compliance can be inherent in 
tubular coolant passage construction, because this contour allows relatively free expansion. Free 
expansion accommodates and distributes the strains preventing local accumulation, and therefore 
extends the life of the chamber. This initial investigation is a simplified, relative comparison of 
the expected lives of milled cooling channels and tubular cooling passages. The analysis compares 
only plug-nozzle chambers that have been used at NASA Lewis Research Center for thrust chamber study 
and development (Fig. 1, Ref. 1). The plug-nozzle chamber has been tested and studied since the 
early 1970’s and the milled channel configuration has been very well characterized. This investiga- 
tion compares two plug-nozzle chamber configurations: (l) a milled OFHC copper channel type 
(Fig. 2); and (2) an OFHC copper tube bundle type (Fig. 3). The comparison is done at the same pres- 
sure loading (i.e., the same chamber pressure, and same coolant pressure) and the same thermal load- 
ing so that the pure compliance, or geometric benefit can be determined. 



FIGURE 2. - CHANNEL PLUG-NOZZLE CHAMBER. 
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FIGURE 3. - TUBULAR PLUG-NOZZLE CHAMBER. 


ANALYTICAL PROCEDURE 

The analytical procedure consisted of a thermal analysis, a strain analysis, and a life correla- 
tion. This three part analysis was performed on a section through the throat of the thrust chamber. 
This area was selected for the analysis because the area of highest heat flux generally produces 
the highest thermal strains. The thermal analysis was performed using SINDA '85, a general purpose 
heat transfer/fluid program that simulates transient or steady-state energy storage and flow in a 
system modeled as a "thermal circuit." The steady state temperature profiles were determined for 
both the channel and tube structures. The channel model used 185 elements, and the tube model used 
180 elements (Figs. 4(a) and (b)). The coolant temperature in both models was set at -360 °F to 
account for the warming of the hydrogen before it gets to the throat. The heat flux for the channel 
model was set at 33 Btu/in.“/sec (Ref. 1). The maximum heat flux for the tube model was 33 Btu/ 
in. 2 /sec and decreased according to the profile reported in Ref. 2. This compensates for the thick- 
ening boundary layer about the periphery of the tube. The coolant-side film coefficient was varied 
in both models until a maximum hot-gas-side wall temperature of 1000 °F was obtained. (Ref. 1 shows 
that 1000 °F was the hot-gas-side wall temperature for past testing.) 

The strain analysis was performed using a two-dimensional elastic-plastic MARC analysis. MARC 
is a general purpose finite element analysis program. Both strain models contained the same number 
of elements as their thermal models. The models consisted of 8-node, three-dimensional, isoparame- 
tric elements (MARC element 7). The strain analysis used the temperatures from the SINDA '85 
analysis along with a chamber pressure of 600 psi and a coolant pressure of 1200 psi . (These are 
standard operating conditions for plug-nozzle tests.) (Ref. 1) The constraints were the same for 
both configurations. The cool outer structure was fixed and the nodes along the planes of symmetry 
assigned zero displacement perpendicular to these planes, but permitted to move in the radial direc- 
tion (Figs. 4(a) and ( b ) ) . The material properties such as coefficient of thermal expansion 
(Fig. 5(a), Ref. 3), yield strength, and modulus of elasticity were input as a function of tempera- 
ture (Figs. 5(b) and (c), Ref. 4). The strain analysis was done using two simplifying assumptions 
and is considered valid when making only relative comparisons. First, no creep analysis was per- 
formed. This absence of creep is considered a valid assumption for the analysis of the plug-nozzle 
chambers because of the short cycle duration (3.5 sec) (Ref. 1). Second, the first cycle was taken 
to be representative of all cycles to follow. The critical strains used for the life analysis were 
determined from taking into account the magnitude of the strains (from the finite element analysis) 
along with other effects such as "ratcheting" (from test experience) that are known to shorten life 
(Ref. 5). These critical strains at steady state conditions were used to determine life instead of 
the more complicated, transiently analyzed strain ranges used in other predictions. 

The life correlation was done using experimental low cycle fatigue data (Fig. 6, Ref. 6). This 
portion of the analysis was not meant to give an absolute value of life, rather it uses the low cycle 
fatigue data as a measure of the material's sensitivity to strain, knowing the materials sensitivity 
to strain, and the relative decrease in strain, a life improvement prediction can be made (Refs. 7 
and 8) . 
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FIGURE 4. - MODELS WITH CONSTRAINTS. 




(C) MODULUS OF ELASTICITY VERSUS TEMPERATURE. 
FIGURE 5. - CONCLUDED. 



ATURE. 

FIGURE 5. - OFHC COPPER MATERIAL PROPERTIES. 



CYCLES TO FRACTURE 


FIGURE 6. - TYPICAL LOW-CYCLE FATIGUE OF OFHC COPPER 
ANNEALED CONDITION AT £ = 0.002 SEC -1 IN ARGON AT 
1000 °F. 
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RESULTS 


This analysis of the two chamber coolant passage configurations was undertaken to compare their 
cyclic lives. This effort began with the thermal analysis using S1NDA '85. The SINDA '85 analysis 
resulted in steady state temperature profiles that are shown in Figs. 7 and 8. The maximum hot- 
gas-side wall temperature was set to 1000 °F and at this condition the profiles were very similar. 

The hottest elements were located at the center of the channel and at the crown of the tube. 

With the thermal profiles, the material properties (as a function of temperature) and test pres- 
sures the strain analysis was performed. This resulted in the strain profiles of the two configura- 
tions as shown in Figs. 9 and 10. These strains are represented as Von Mises strains which are 
equivalent uniaxial strains. As can be seen from the strain contour (Fig. 9) for the channel struc- 
ture, the location of the highest strain is at the center of the rib. This is not the typical fail- 
ure site for this structure when chambers are actually tested, Fig. 11 shows an actual failure that 
has the rupture at the center of the cooling passage. This is not considered to be in error because 
in actual testing there are ratcheting affects (Ref. 9) taking place at the center of the channel 
which reduce the life dramatically (Ref. 5). The strain at the center of the channel is just 
slightly less than the maximum strain at the rib, but when combined with the racheting effect the 
critical strain area is at the center of the channel. 


The strain contour (Fig. 10) for the tube structure shows that the highest strain occurs at the 
crown of the tube. The failure is expected to occur at the center line of the tube crown for the 
same reason as the channel above, high strain and racheting. 


Figure 12 is a summary of the life analysis performed. The plot indicates experimental data of 
life as a function of hot-gas-side wall temperature for milled 0FHC copper channel chambers (Ref. 1). 
It also shows the calculated life versus hot-gas-side wall temperature for the same configuration. 
There is very good agreement between the experimental data and the calculated data. The calculated 
data is slightly higher because of the simplifying assumptions made for the analysis and the fact 
that the number of cycles to failure was determined using low cycle fatigue specimen data. The low 
cycle fatigue specimen data is the average life at a particular strain range. The low cycle fatigue 
curve is based upon a series of data points which have a degree of scatter, consequently the line 
represents some average. Experimentally the plug-nozzle chambers are considered failed when the 
first sign of a crack through the wall occurs. These chambers have 72 channels, hence 72 fatigue 
specimens. Therefore, when the first crack occurs through the wall the weakest specimen is taken as 
the data point not the average as in the low cycle fatigue data. Figure 12 also shows the calcula- 
ted life of the tube as a function of hot-gas-side wall temperature, showing that the tube cooling 
passage has promise of living approximately 100 percent longer than the channel cooling passage. 
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FIGURE 7. - CHANNEL MODEL 
TEMPERATURE CONTOURS. 


FIGURE 8. - TUBE MODEL 
TEMPERATURE CONTOURS. 


FIGURE 9. - CHANNEL MODEL FIGURE 10. - TUBE MODEL 

STRAIN CONTOURS. STRAIN CONTOURS. 
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FIGURE 11. - TYPICAL CHANNEL WALL FAILURE. 



FIGURE 12. - HOT-GAS-SIDE WALL TEMPERATURE VERSUS 
CYCLES TO FAILURE. 


SUMMARY OF RESULTS 

An increase in life of a tubular coolant passage configuration, when compared to an equivalent 
channel configuration, of ~100 percent was found for hot-gas-side wall temperatures ranging from 
800 to 1000 °F . A summary of critical strains along with the corresponding stresses and number of 
cycles to failure are provided in Fig. 13. A comparison of the number of cycles to failure as a 
function of hot-gas-side wall temperature for the channel and tube cooling passages is also provided 
in Fig. 13. 



HOT-GAS- 
SIDE WALL 
TEMPERATURE, 
°F 

STRESS, 

PSI 

CRITICAL 
STRAIN, 
IN. /IN. 

CYCLES TO 
FAILURE 

MILLED 

CHANNEL 

PASSAGE 

1000 

8850 

0.0125 

350 

900 

8025 

0.0108 

500 

800 

8001 

0.0091 

740 

TUBE 

PASSAGE 

1000 

6969 

0.0094 

700 

900 

7081 

0.0083 

900 

800 

7184 

0.0073 

1200 


FIGURE 13. - SUMMARY TABLE OF HOT-GAS-SIDE WALL TEMPERATURE, 
STRESS, STRAIN, AND CYCLES TO FAILURE. 


RECOMMENDATIONS 

Three recommendations are made for tubular configuration life analysis. First, although this 
analysis is valid for relative comparisons of channel and tube structures, a three-dimensional visco- 
plastic analysis of the tube structure is required to get a accurate understanding of what is happen- 
ing. Second, alloy tubes should be analyzed to determine if alloy properties are applicable when 
compliance is the goal. Third, the results of both this analysis and subsequent analyses need to be 
validated with testing. 
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